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1 Introduction 

Mathematics is a summary and abstraction of the production practices and the natural 
sciences, and is a powerful tool to explain natural phenomena and reveal the laws of nature 
as well. Atmospheric circulation is one of the main factors affecting the global climate, so it 
is very necessary to understand and master its mysteries and laws. Atmospheric circulation 
is an important mechanism to complete the transports and balance of atmospheric heat 
and moisture and the conversion between various energies. On the contrary, it is also the 
important result of these physical transports, balance and conversion. Thus it is of necessity 
to study the characteristics, formation, preservation, change and effects of the atmospheric 
‘Foundation item: the National Natural Science Foundation of China (NO. 11271271). 
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circulation and master its evolution law, which is not only the essential part of human’s 
understanding of nature, but also the helpful method of changing and improving the accuracy 
of weather forecasts, exploring global climate change, and making effective use of climate 
resources. 

The atmosphere and ocean around the earth are rotating geophysical fluids, which are also 
two important components of the climate system. The phenomena of the atmosphere and 
ocean are extremely rich in their organization and complexity, and a lot of them cannot be 
produced by laboratory experiments. The atmosphere or the ocean or the coupled atmosphere 
and ocean can be viewed as an initial and boundary value problem [8] [10], or an infinite 
dimensional dynamical system [3] [4] [5]. These phenomena involve a broad range of temporal 
and spatial scales [2]. For example, according to [11], the motion of the atmosphere can be 
divided into three categories depending on the time scale of the prediction. They are motions 
corresponding respectively to the short time, medium range and long term behavior of the 
atmosphere. The understanding of these complicated and scientific issues necessitate a joint 
effort of scientists in many fields. Also, as [11] pointed out, this difficult problem involves a 
combination of modeling, mathematical theory and scientific computing. 

Some authors have studied the atmospheric motions viewed as an infinite dimensional 
dynamical system. In [3], the authors study as a first step towards this long-range project 
which is widely considered as the basic equations of atmospheric dynamics in meteorology, 
namely the primitive equations of the atmosphere. The mathematical formulation and at¬ 
tractors of the primitive equations, with or without vertical viscosity, are studied. First of 
all, by integrating the diagnostic equations they present a mathematical setting, and obtain 
the existence and time analyticity of solutions to the equations. They then establish some 
physically relevant estimates for the Hausdorff and fractal dimensions of the attractors of 
the problems. In [12], based on the complete dynamical equations of the moist atmospheric 
motion, the qualitative theory of nonlinear atmosphere with dissipation and external forcing 
and its applications are systematically discussed by new theories and methods on the infinite 
dimensional dynamical system. In [13], by Lions theorem in the Hilbert space, the existence 
and uniqueness of the weak solution of water vapour-equation with the first boundary-value 
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problem are proven, and the scheme of the finite-element method according to the weak 
solution is proposed. In [14], the model of climate for weather forecast is studied, and the 
existence of the weak solution is proved by Galerkin method. The asymptotic behaviors of 
the weak solution is described by the trajectory attractors. 

In this article, Atmospheric circulation equations with humidity effect is considered, which 
is different from the previous research. Previous studies are based on two kinds of equations, 
one is about the heat and humidity transfer([12],[13]), without considering the diffusion of 
heat and of humidity; the other is p-coordinates equation in([3], [14]), being used to consider 
the horizontal movement of the atmosphere, and the vertical direction of the velocity is 
transformed into pressure and topography. Atmospheric circulation equations with humidity 
effect, derivation from the original equations in [9], considering the diffusion of heat and of 
humidity, not be restricted by p-coordinates, can be deformed based on the different concerns. 
In the last part of the article, Existence of global solutions to the atmospheric circulation 
models is obtained, which implies that atmospheric circulation has its own running way 
as humidity source and heat source change, and confirms that the atmospheric circulation 
models are reasonable. 

The paper is organized as follows. In Section 2 we present derivation of atmospheric 
circulation models. In Section 3, we prove that the atmospheric circulation models possess 
global solutions by using space sequence method. 

2 Derivation of Atmospheric Circulation Models with 
Humidity Effect 

2.1 Original Model 

The hydrodynamical equations governing the atmospheric circulation are the Navier-Stokes 
equations with the Coriolis force generated by the earths rotation, coupled with the first law 
of thermodynamics. 

Let (ip, 9, r) be the spheric coordinates, where <p represents the longitude, 6 the latitude, 
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and r the radial coordinate. The unknown functions include the velocity field u = (u^, ug, u r ), 
the temperature function T, the humidity function q, the pressure p and the density function 
p. Then the equations governing the motion and states of the atmosphere consist of the mo¬ 
mentum equation, the continuity equation, the first law of thermodynamic, and the diffusion 
equation for humidity, and the equation of state (for ideal gas): 

p[^ + + 2 ft x u] + Vp + itpg = pAip (2.1) 


+ div(pu) = 0, 


pcu[— + u ■ VT] + pdivu = Q + ktAT, 


Pi ^ + u ■ Vg] = G + u q Aq, 


P = RopT, 


Where —oo < p < +oo, — ^ < 9 < ro < r < ro + h, ro is the radius of the earth, h is 
the height of the troposphere, Q is the earth’s rotating angular velocity, g is the gravitative 
constant, p, kt, K q , c v , Rq are constants, Q and G are heat and humidity sources, and 1 1 = 
(0, 0,1). The differential operators used are as follows: 

1. The gradient and divergence operators are given by: 


V ^ d 1 d d 
r cos 9 dp ’ r 89' dr 

l9.o. 1 d(ug cos 9) 

i,VU = ^Wr {r Ur) + r cos 6 80 


1 dug, 

r cos 9 dp 


2. In the spherical geometry, although the Laplacian for a scalar is different from the 
Laplacian for a vectorial function, we use the same notation A for both of them: 


. . . 2 du r 2 sin 9 dug 

Au = (A Up ~\ 2 --1 2 - To~a - 2 - 27 )’ 

r z cos 9 dp r z cos z 9 op r z cos z 9 


Aug + 


A u r + -AL + 


2 du r ug 2sin# Up 

r 2 89 r 2 cos 2 9 r 2 cos 2 6 dp ’ 

2 d(ugcos 9) 2 Up. 

r 2 r 2 cos 9 89 r 2 cos 9 dp ’ 
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A f _ 1 df q<L\ a 1 d 2 / 1 df 2 9 

r 2 cos 9 89 89 r 2 cos 2 0 dtp 2 r 2 3r <9r 

3. The convection terms are given by 


V7 / v-7 , UpU r U^Ug 

V u u = (u ■ Vu v H— 21 - 1 — tan 0, 


u ■ vug H-1— 11 tan#, u ■ vu r - r - 


4. The Coriolis term 2 ft x u is given by 


2 ft x u = 2fl(cos 0u r — sin Oug, sin Oug , — cos 0-u 


<PJ> 


Here ft is the angular velocity vector of the earth, and f l is the magnitude of the angular 
velocity. 

They are supplemented with the following initial value conditions 


(u,T,q) = (ip w ,<p2o,V3o) at t = 0. 


( 2 . 6 ) 


Boundary conditions are needed at the top and bottom boundaries (ro,ro + h). At the 
top and bottom boundaries (r = ro, ro + h), either the Dirichlet boundary condition or the 
free boundary condition is given 


(Dirichlet) 


(u,T,q) = (0,To,go)i r = r 0 , 
(u,T,q) = (0,Ti,qi), r = r 0 + h, 


(2.7) 


(free) 


(u,T,q) = (0 ,T 0 ,q 0 ), 


r = r 0 , 


(ur,T,q) = ( 0,Ti,gi), a (yo) = q r = rQ + h, 
For p, periodic condition are usually used, for any integers k\ € Z 

(u,T,q)(ip + 2kiir,0,r) = (u,T, q)(<p + 2kiir,Q, r). 


( 2 . 8 ) 


(2.9) 


Because (<p, 0, r ) are in a circular field with C°° boundary, the space domain is taken as 
(0, 27 t) x (—f, f) x (ro, ro + h ) and periodic condition is written as 


(u, T, g)(0, 0, r) = (u, T, q)(2n, 0, r ). 


( 2 . 10 ) 
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For simplicity, we study the problem with the Dirichlet boundary conditions, and all re¬ 
sults hold true as well for other combinations of boundary conditions. Atmospheric convection 
equations can be read as (2.1)-(2.7), (2.10). 

The above equations were basically the equations used by L. F. Richardson in his pioneer¬ 
ing work [9]. However, they are in general too complicated to conduct theoretical analysis. 
As practiced by the earlier workers such as J. Charney, and from the lessons learned by 
the failure of Richardson’s pioneering work, one tries to be satisfied with simplified models 
approximating the actual motions to a greater or lesser degree instead of attempting to deal 
with the atmosphere in all its complexity. By starting with models incorporating only what 
are thought to be the most important of atmospheric influences, and by gradually bring¬ 
ing in others, one is able to proceed inductively and thereby to avoid the pitfalls inevitably 
encountered when a great many poorly understood factors are introduced all at once. The 
simplifications are usually done by taking into consideration of some main characterizations 
of the large-scale atmosphere. One such characterization is the small aspect ratio between the 
vertical and horizontal scales, leading to hydrostatic equation replacing the vertical momen¬ 
tum equation. The resulting system of equation are called the primitive equations (PEs); see 
among others [3]. The another characterization of the large scale motion is the fast rotation 
of the earth, leading to the celebrated quasi-geostrophic equations [1]. 

For convenience of research, the approximations we adopt involves the following compo¬ 
nents: 

First, we often use Boussinesq assumption, where the density is treated as a constant 
except in the buoyancy term and in the equation of state. 

Second, because the air is generally compressible, we do not use the equation of state 
for ideal gas, rather, we use the following empirical formula, which can by considered as the 
linear approximation of 

P = Po[l - «t(P - T 0 ) + a q (q - g 0 )], (2.11) 

where po is the density at T = To and q = qo, and ax and a q are the coefficients of thermal 
and humidity expansion. 
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Third, since the aspect ratio between the vertical scale and the horizontal scale is small, 
the spheric shell the air occupies is treated as a product space S% 0 x (ro, ro + h). This approxi¬ 
mation is extensively adopted in geophysical fluid dynamics. Under the above simplification, 
we have the following equations governing the motion and states of large scale atmospheric 
circulations: 

du \ 1 

— + V„tt = vAu -2 Q xu - Wp - [1 - a T (T - T 0 ) + a q (q - q 0 )\ge z , (2.12) 

dt p 0 


+ ( u • V)T — Q + ktAT, 
^ + (u ■ V)q = G + n q Aq, 


(2.13) 

(2.14) 


divu = 0 , 


(2.15) 


where (p, 0 , z) £ M = S 2 0 x (ro, ro + h), v = A i s the kinematic viscosity, u = u v e v + ugeg + 
u r e r , (e^,e 0 ,e r ) the local normal basis in the sphereric coordinates, and 


V u it = ((it • V)u„ -)— 1 -^-^tan6*)e v + 

ro r 0 


uqu z it; 


((it • V)u <9 H- 21 + — tan 0)e<j + ((it • V)it z - 


ro ro 


u% + u 2 e 
ro 


) e z, 


. , . 2 9tt z 2 sin 6 dua Um 

Alt — (A U p + j 2 _ n ""o "7 + __ 2 _ 9 . /i T 7 _2_ 9 /i) e ¥> + 


cos 0 dip Tq cos 2 0 dp 7q cos 2 0 


. . 2 9it z 

(Aitg -|—2 - 


ue 


2 sin 0 it, 


t-q 90 rg cos 2 0 cos 2 0 9</? 


)ee+ 


, A-i- | 2u 0 , 2 9(it0 cos 0) 2 

\2±Uz T" o > . 2_/i o/i V_/i o _/®2> 


r 2 Tq cos 0 90 


Tq cos 0 dp' 

l dp 1 dp dp 

vp = -+ jr e zi 

ro cos 0 90 r 90 9z 

1 du P 1 9(it0cos0) du z 

ro cos 0 9</? ro cos 0 90 dz ’ 

and the differential operators (it • V) and A are expressed as 

f V7N % 9 | Ug d t d 

(it • V) — --r 7: -1-— + U z J— , 

ro cos 0 dp ro 90 9z 

a 1 <9 2 1 9 9 9 2 

A = —-+ o--Wx(cOS0 —) + T£J2- 


rg cos 0 9 </? 2 r§ cos 0 90 v dO dz 2 
Equations (2.12)-(2.15) are supplemented with boundary conditions (2.7), (2.10). 
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2.2 Simplification of Model 


Atmospheric circulation is the large-scale motion of the air, which is essentially a thermal 
convection process caused by the temperature and humidity difference between the earths 
surface and the tropopause. It is a crucial means by which heat and humidity are distributed 
on the surface of the earth. Air circulates within the troposphere, limit vertically by the 
tropopause at about 8-lOkrn. Atmospheric motion in the troposphere, together with the 
oceanic circulation, plays a crucial role in leading to the global climate changes and evolution 
on the earth. There are two types of circulation cells: the latitudinal circulation and the lon¬ 
gitudinal circulation. The latitudinal circulation is characterized by the Polar cell, the Ferrel 
cell, and the Hadley cell, which are major players in global heat and humidity transport, and 
do not act alone. The zonal circulation consists of six circulation cells over the equator. The 
overall atmospheric motion is known as the zonal overturning circulation, and also called 
the Walker circulation. The most remarkable feature of the global atmospheric circulation 
is that the equatorial Walker circulation divides the whole earth into three invariant regions 
of atmospheric flow: the northern hemisphere, the southern hemisphere, and the equatorial 
zone. We also note the important fact that the large-scale structure of the zonal overturning 
circulation varies from year to year, but its basic structure remains fairly stability, it never 
vanishes. Based on these natural phenomena, we here present the Zone Hypotheses for at¬ 
mospheric dynamics, which amounts to saying that the global atmospheric system can be 
divided into three sub-systems: the North-Hemispheric System, the South-Hemispheric Sys¬ 
tem, and the Tropical Zone System, which are relatively independent, and have less influence 
on each other in their basic structure. More precisely, the Atmospheric Zone Hypothesis is 
stated in the following form[6]. 

Atmospheric Zone Hypothesis. The atmospheric circulation has three invariant re¬ 
gions: the northern hemisphere domain (0 < 6 < ^), the southern hemisphere domain 
(—^ < 0 < 0), and the equatorial zone (0 = 0). Namely, the large-scale circulations in their 
invariant regions can act alone with less influence on the others. In particular, the velocity 
field u = (u^.ug^Uz) of the atmospheric circulation has a vanished latitudinal component in 



a narrow equatorial zone, i.e., ug = 0, for —e < 9 < e, where e > 0 is a small number. 

Atmospheric Zone Hypothesis is based on the following several evidences from theory and 
practice: 

(1) The global atmospheric motion equations (2.12)-(2.15) are of 9— reflexive symmetry, 
i.e. under the 9— reflexive transformation (p,9,z) —> ((p,—9,z), the velocity field u becomes 
(u v ,ug,u z ) —y (up, — ug, u z ), and equations (2.12)-(2.15) are invariant, which implies that 
this system is compatible with the Atmospheric Zone Hypothesis. 

(2) Climatic observation data show that when the El Niho-Southern Oscillation (the 
behavior that the Walker circulation cell in the Western Pacific stops or reverses its direc¬ 
tion) takes place, no oscillation occurs in the latitudinal cells. It demonstrates the relative 
independence of these circulations in their invariant domain. 

(3) When a cold current moves southward from Siberia, or a violent typhoon sweeps north¬ 
ward from the tropics, the weather in Southern Hemisphere has no response. Atmospheric 
Zone Hypothesis provides a theoretic basis for the study of atmospheric dynamics, by which 
we can establish locally simplified models to treat many difficult problems in atmospheric 
science. 

(4) For the three-dimensional atmospheric circulation equation, it is too difficult to study. 
So we study the equatorial atmospheric circulation. 

The atmospheric motion equations over the tropics are the equations restricted on 9 = 0, 
where the meridional velocity component ug is set to zero, and the effect of the turbulent 
friction is taking into considering 


. U^nUz 

— (u • VK - - - 1 - v(Autp + 


rl dp 


2U V N OO 1 d P 

—J 1 ) ~ ootv - 2 Sluz - 

rt Po r o dip 


(2.16) 


— 1 = -{u ■ V)u z + — + u(Au z + - 2 ~) - U z - 2 Qup 

dt ro rg dip rg 


1 dp 
Po dz 


[1 - a T {T - T 0 ) + a q {q - q 0 )\g, 


(2.17) 


dT 

~dt 


— (u • V)T + /vpAT + Q , 


(2.18) 
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— = -(u ■ V)<? + n q Aq + G, 


(2.19) 


1 duy du^ _ 

ro dip dz ’ 


( 2 . 20 ) 


Here <7* = Cih 2 {i = 0,1) represent the turbulent friction, ro is the radius of the earth, the 
space domain is taken as M = x (ro,ro + h) with S} Q being the one-dimensional circle 


with radius ro, and 


, u v d d . id 2 d 2 

( “' V)_ ^ + “ 2 d? A ~7 2 W I+ d^ 


For simplicity, we denote 


(xi,x 2 ) = ( r 0 <p,z ), (ui,u 2 ) = (u<p,u z ). 


The atmospheric motion equations(2.16)-(2.20) can be written as 


dm . . «i« 2 . / A ,2 du 2 2«i on 1 d P 

= ~(u ■ v)u\ -h v(Aui H---o") - ctqU\ - 2ilu 2 -w—, 

ot ro ro oxi po ox\ 

du 2 , iii , . 2 du\ 2u 2 . ^ 

= —(it • V)u 2 + — + v{Au 2 + ——- 2 ") - criit 2 - 2Q,u\ 

ot ro ro ox\ r q 


( 2 . 21 ) 


1 dv 

~Vo~dx~2 ~ ^ ~ ~~ + aq ^ ~ % ^ 9 ' 


( 2 . 22 ) 


= -( m .V)T + « t AT + Q, 
dt 


(2.23) 


— = -(u- S7) q + Kq Aq + G, 


(2.24) 


dui du2 _ 
dx i dx 2 


(2.25) 


To make the nondimensional form, let 


, / h 2 i kt ' 

x = hx , t = —t u = —it , 
kt h 


T = T 0 - (T 0 - Tx)-^ + (T 0 - T X )T , 

q = qo - (qo - n)y + (qo - qi)q', 
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PQVKTP r . «T, t «<Z/ x*2i 

P= —^2 -+ ^-(To-Ti)--—(go - QiJ-f^U 

The nondimensional form of (2.21)-(2.25) reads 

du\ i i v A / cr 0 ^ 2 ' 2f lh 2 / 1 / Op 7 

= _( u . V)«! + —Auj-«!-w 2 - T^-r 

ot Kj 1 K>t Kt /c T OXI 


h 2 h du 2 2 h? 

— u\u 2 H- 

n> 


rokr dx\ ro kt^ 1 ’ 


(2.26) 


3^2 

dt' 


, i . / i/ / cri/i 2 / 2fl/r 2 / h 
-(u • \/)u 2 H-Au 2 -— 


2 h du, 2hr / 
Mi-Ml H-— —;-O Uo 


Kt 


- “2 -“i “i i-; ttu -o— 

act act ro ro^T o,t 2 a-qACt 


+ ^-[a T (T° - Ti)T - ag(g° - gi)g], 


act 3^2 AC 7 


(2.27) 


<9T 


/r 2 


8f r “ ' V)r + AT + “ 2 + (To -TOkt 0, 


(2.28) 


3a/ . / / Ac a . / / /i 2 

_ = . V) , + -A„ +-2 + (To _ TiKt 


G, 


(2.29) 


divu = 0. 


(2.30) 


Because ro is far larger than iti, u 2 , the atmospheric motion equations (2.26)-(2.30) can 
be read 


3't/, 

3^ 


——ttt + ^-[olt{Tq - Ti)T' - a q (q 0 - qi)q], (2.32) 

act ox 2 ACy 

37~’ / /? 2 

lr = -(„'.V)T' + AT' + „; + M -^<J, (2.33) 

W = “ ( “ ' V)9 ’ + S A? ' + “ 2 + (T„ - T,) W G ’ <2 ' 34) 


/ / V k ' °- 0 /l 2 / 212/l 2 / 1 / 3p 7 

-(u • V) , u 1 H- Amj- u -2 - u 2 - —r, 

kt kt kt act ox 1 


(2.31) 


du 2 , / . / !A / (Ti/a 2 / 

—^ = -(u • V)u 2 H -Am 2 -Wo - 

ot kt kt 


2Qh 2 

kt 


-Wi 
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divu = 0 . 


( 2 . 35 ) 


Let P r — Kt , L e 

Q' = I= 

(2.35) reads 


_ K g o gar (T°—T\)h? 5 ga q (q°-qi)h 3 ' <7.;fe 2 2Qh 2 

K T > 21 > 21 K'J'V ' u i V 1 W V 1 

(Tq—Ti)k. t C om itting the primes, the nondimensional form of (2.31)- 


P r (Au — Xp — au) + P r (RT — Rq)iZ — (V • u)u, 

(2.36) 

dT 

= AT + u 2 — (u ■ V)T + Q , 

(2.37) 

da 

— = L e Aq T u 2 (u X)q + G, 

(2.38) 

divu = 0, 

(2.39) 


where a is constant matrix 

/ cr 0 u 

a = 

\ UJ <7l 

The problem (2.36)-(2.39) are supplemented with the following Dirichlet boundary con¬ 
dition at x 2 = 0,1 and periodic condition for x\\ 


(u,T,q) = 0, x 2 = 0,l, 


(2.40) 


(u,T,q)( 0,x 2 ) = (u,T,q)(27r,x 2 ), 


(2.41) 


and initial value conditions 


(u,T,q) = (u 0 ,T 0 ,qo), t = 0. 


(2.42) 


3 Global Solution of Atmospheric Circulation Equations with 
Humidity Effect 

3.1 Preliminaries 

Let X be a linear space, X\, X 2 two separable reflexive Banach spaces, H a Hilbert space. 
Xi, X 2 and H are completion space of X under the respective norms. X\, X 2 C H are 
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dense embedding. F : X 2 x (0, oo) —> A* is a continuous mapping. We consider the abstract 
equation 


= Fu, 0 < t < oo, 


[ u(0) = tp, 

where ip £ H, u : [0, +oo) — > FI is unknown. 

Definition 3.1 Let ip £ H be a given initial value, u £ L p ((0, T), X 2 ) D L°°((0,T), ff), 
(0 < T < 00 ) is called a global solution of Eq(3.1), if u satisfies 

(u(t),v)n = f < Fu,v > dt + (ip,v)n, Vu € X\ c H. 

Jo 

Definition 3.2 Let u n , uq £ L p (( 0, T),X 2 ). u n uq is called uniformly weak convergence 
in L p ((0,T), X 2 ), if {u n } C L°°((0 ,T),H) is bounded, and 


U n 1 


in I7((0,T),A 2 ), 


{ lim, woo / 0 T | < u n - uq, v > H 1 2 dt, = 0, \/v £ H. 

Definition 3.3 A mapping F : X 2 x (0, 00 ) —> AJ is called T-weakly continuous, if for 
p = (pi ,..., p m ), 0 < T < 00 , and u n uniformly weakly converge to uq under Eq(3.2), we 


lim [ <Fu n ,v>dt= I <Fuo,v>dt, Mv £ X\. 

n ^°° Jo Jo 

Lemma 3.4 Let {n n } C T p ((0, T), W m,p )(m > 1) be bounded sequence, and {rt n } uni¬ 
formly weakly converge to uo C L p ((0,T) 1 W m,p ) 1 i.e. 

f u n -+u 0 in L p ((0,T),W m ' p ), p> 2, 


[ lirn^oo / 0 T [J n (u n - u 0 )vdx] 2 dt = 0, Vv £ 
Then for all |a| < m — 1, we have 


D a u n ^D a u 0 in L 2 ((0, T) x 12). 


Lemma 3.5 ^ Assume F : X 2 x (0, 00 ) —> X\ is T-weakly continuous, and satisfies 
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(Al) there exists p = (pi,- ■ ■ ,p m ), Pi > 1(1 < i < m), such that 

< >< — Ci||it||j ^ 2 + C2 + f(t), Vw € A, 

where C\,C 2 are constants, / € L 1 (0,T)(0 <T < 00 ), || • = Y^IL\ | • |f% | * |j is seminorm 

of A 2 , || • ||.y 2 = E™i I • U, 

(A2) there exists 0 < a < 1, for all 0 < h < 1 and u £ C' 1 ([0, 00 ), X), 

rt+h 

| j < Fu,v > dt | < Ch a , Mv £ X and 0 <t <T, 

where C > 0 depends only on T, ||u||xi, fg IMIx 2 ^ an( l su P 0 <t<T IMIff- 
Then for all (p £ H, Eq(3.1) has a global weak solution 

u £ L°°((0,T),H) nL p ((0,T),X 2 ), 0 < T < 00 , p in (Al). 

Remark 3.6 || • ||x denotes norm of X, and C\ are variable constants. 

3.2 Existence of Global Solution 

We introduce spatial sequences 

X = {</>= (u,T,q) £C°°(n,R 4 )\(u,T,q) satisfy (2.39) - (2.41)}, 

H = {(/) = (u,T,q) £ L 2 (Q,R 4 )\(u,T,q) satisfy (2.39) - (2.41)}, 

H x = {(/) = {u, T, q) £H 1 {n,R 4 )\(u 1 T,q) satisfy (2.39) - (2.41)}, 

Theorem 3.7 If cj)g = (uo,To,qo) £ H, Q,G £ L 2 (Q), then Eq(2.36)-(2.42) there exists a 
global solution 

(u,T,q) £ L°°((0, T),H) D T 2 ((0, T),H\), 0 < T < 00 . 

Proof. Definite F : H\ —> H\ as 

< F(j),iJj > = / n [— P r VnVv — P r au ■ v + P r (RT — Rq)v 2 — ( u • V)u • v 

-VTS7S + u 2 S - (u ■ V)TS + QS - L e VqVz + u 2 z 

— (u ■ V)qz + Gz\dx, = (v, S, z) £ H\. 
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Let 'ip = (j). Then 


< F(j>, <f> > = / q [— P r \\7u\ 2 — P r au ■ u + P r (RT — Rq)u 2 — (u ■ V)it • u 

— |VT| 2 + u 2 T - (■ u • V)TT + QT - L e \Vq\ 2 

+u 2 q — (u ■ V)qq + Gq]dx 

= / q [— P r \Vu\ 2 — P r au ■ u + ( P r R + 1)u 2 T — (. P r R — l)qu 2 

— |VT| 2 + QT- L e \Vq\ 2 + Gq]dx 

< -Ci / n [|Vu| 2 + |VT| 2 + \Vq\ 2 \dx + C 2 / n [|u| 2 + |«||T| + |g||u 


+ \Q\\T\ + \G\\q\}dx 

< -Cl JJVul 2 + |VT | 2 + \Vq\ 2 ]dx + C 2 / n [|n | 2 + |T | 2 + \q\ 2 ]dx 
+ ( ^ 3 /n[|Q| 2 + \G\ 2 ]dx 

< —Ci\\(p\\ 2 Hl +C 2 \\cp\\ 2 H + Ci, 

which implies (j4i). 

For V</> € L 2 (0,T), Hi) n L°°((0, T),H) and ip£X,h{0<h< 1), we have 
I It +h < F<f>,ip > dt\ 

= I It t+h In[-PrVuVv - P r au • v + P r (RT — Rq)v 2 — (u • V)tt • v - VTXS 


+u 2 S — (u ■ V)TS + QS — L e X7qVz + u 2 z — (u ■ \7)qz + Gz\dxdt\ 
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< C'/ t m /n[|V«||V^| + |u||v| + \T\\v\ + |?||t;| + |VT||VS| + M|S| 

+ |<2||S'| + |Vg||Vz| + [u\\z\ + \G\\z\]dxdt + C // +ft [| f n (u ■ V)u ■ vdx\ 

+ | Jq(u ■ X7)TSdx\ + | fft(u ■ \7)qzdx\]dt 

< CJt +h [(Jn \Vu\ 2 dx)^(J n |Vu| 2 dx)^ + (f a \T\ 2 dx)^{f n \v\ 2 dx)^ 

+(/n \q\ 2dx )Hln \v\ 2 dx)^ + (f n \VT\ 2 dx)^{J n \VS\ 2 dx )3 
+(/n l«l 2 *«0*(/n \S\ 2 dx)t + dn \Q\ 2 dx)Hh I S\ 2 dx)h 

+(/n |V?| 2 dx)^(J n \Vz\ 2 dx )s + (f Q \u[ 2 dx)5(f n \z\ 2 dx)^ 

+(S {2 \ G \ 2 dx)Hln \z\ 2 dx)^]dt + C J t t+h [\ Ei,j=i Jnhuj^-dx\ 

+| Eh In UiT^dx | + | Eh fn w£dx\]dt 

< C , (II' u IIl 2 (0,T;_H' 1 )II^ ) ' 1; IIl 2 ^ 2 + I|-^'IIl 2 (0,T;L 2 ) II^IIl 2 ^' 2 + ll^/ll L 2 (0,T;L 2 ) 11^IIL 2 ^ 2 

+ II^IIl 2 (0,T;H 1 )II^£||l 2 ^ 2 + IMIl^(0,T;L 2 )I|£||l 2 ^ 2 + IIQIIi 2 II'S’IIl 2 ^' 

+ IMIl 2 ( 0 ,T;H 1 )II^~IIl 2 ^ 2 + \\ u \\L°°( 0 ,T-,H)\\z\\L 2 h + II^11L 2 II 2 IIL 2 h 

1 1 

+ ll' u llc 1 II ^IIl°°( 0,T,L 2 )^' + ll'S'llc' 1 \\u\\ loo (0,T;i 2 ) mi L°° (0,T;L 2 ) ^ 

1 1 

+ lkllc 1 IMI L°°(0,T;L 2 ) IMIl°°(0,T;L 2 )^) 

< Chk, 
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which implies (A 2 ). 

We will prove that F : H\ —»• HI is T-weakly continuous. Let <p n = (■ u n ,T n ,q n ) —^ 
( u°,T°,q °) be uniformly convergence, i.e., {4> n } C L°°((0,T); H) is bounded, and 

I 4>n^4>o in LP((0,T);F!), 

[ lim, woo / 0 T I < 4>n - 00, ip >H 1 2 dt = 0, € H. 


From Lemma 3.4, we known that 4> n —> <fi 0 in L 2 (i7 x (0, T)). 

Then Vip £ X C R 4 ) H H\ . we have 

pt p pt p ^ Oil™ 

lim / / (u n ■ X)u n ■ vdxdt = lim y vF—^-Vjdxdt 

- 1 **-►«> Jo ./q z ^ l * J 


and 


and 


n—too Jq J n 


= — lim 

n—>-oo 


V uyy^p-dxdt 

0 ./otj- 1 3 »■>■> 


— [ [ (u° ■ V)v ■ u°dxdt 

Jo Jn 

[ f (u° ■ V)u° • vdxdt, 

Jo Jn 


n ~^°° Jo Jn 


lim / {u n ■ V)T n Sdxdt = lim Y 


= - lim f 1 f 

n^ooj 0 J n 


j 0 Jn 

OS 
dxi 


,dT" 

dxi 


■Sdxdt 


f)Q 

y vFT n ——dxdt 
2—1 


- T f (u° ■ X7)ST°dxdt 

Jo Jn 

C [ (u° ■ V)T° Sdxdt, 

Jo Jn 


lim 

n—too Jo Jn 


{u n ■ X)q n zdxdt = lim / / ^ u™ —— zdxdt. 


,dq r ‘ 


n ->°° Jo Jn f—J " dxi 


dz 


= - lim / / y v?q n —-dxdt 

L I,, ^ dxi 


o Jn i=1 


- f [ (u° ■ V)zq°dxdt 

Jo Jn 

[ f (u° ■ X)q°zdxdt, 
Jo Jn 
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Thus 


< Fcf> n ,ip > dt = [ <F4>o,^>dt. (3.4) 

Jo 

Because X is dense in Hi, Eq(3.4) holds for ^ G H\. In other words, the mapping 
F : H\ —> HI is T-weakly continuous. 

From Lemma 3.5, Eq(2.36)-(2.42) has a global weak solution 

(u,T,q) € L°°((0, T),H) D L 2 ((0, T),H\), 0 < T < oo. 

□ 

Remark 3.8 Existence of global solutions to the atmospheric circulation models implies 
that atmospheric circulation has its own running way as humidity source and heat source 
change, and confirms that the atmospheric circulation models are reasonable. 
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